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SUMMARY 


The impact motion of the inflated sphere landing vehicle with a 
payload centrally supported from the spherical skin hy numerous cords has 
been determined on the assumption of uniform isentropic gas compression 
during impact. The landing capabilities are determined for a system 
containing suspension cords of constant cross section. 

The effects of deviations in inpact velocity and initial gas 
temperature from the design conditions are studied. Also discussed are 
the effects of errors in the time at which the skin is ruptured. These 
studies indicate how the design parameters should be chosen to insure 
reliability of the landing system. 

Calculations have been made and results are presented for a sphere 
inflated with hydrogen, landing on the moon in the absence of an 
atmosphere. The results are presented for one value of the skin-strength 
parameter . 


INTRODUCTION 


Various devices and techniques are being considered as solutions to 
the problem of alleviating or minimizing the acceleration of a payload 
during impact with the surface of the moon or the planets. Among the 
simpler methods (in concept) which have been proposed is the gas-inflated 
inpact bag. A possible configuration, the inflated sphere landing vehicle, 
was described and its motion and performance were analyzed in references 
1 and 2. 

The inflated sphere landing vehicle has a strong and flexible, but 
relatively nonstretchable skin. The inflating gas is light, but has high 
energy-absorbing properties. The best such gases are hydrogen and helium. 
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The payload package to he landed is suspended in the center of the sphere 
hy numerous cords attached to the skin. Because of the nonstretchahle 
skin and the internal gas pressure, the vehicle maintains a truncated 
spherical shape during the impact as the part of the skin in contact with 
the impact surface is collapsed, folded, and held on the surface hy the 
gas pressure. The payload comes to rest on the surface when the kinetic 
energy of the system has heen absorbed hy the gas in the vehicle. The 
skin may he ruptured at the instant of zero velocity to prevent rehound. 

The results obtained from the analysis of reference 1 showed that the 
maximum impact velocity was limited hy the skin strength to values less 
than about 1000 feet per second. The vehicle could, however, he advanta- 
geously used in the terminal phase of a landing maneuver because it 
functions without the need for maneuverability, guidance, control, attitude 
stabilization, or complex landing gear. 

In the analysis of reference 1 the approximation was made that the 
pressure throughout the inflating gas was uniform at a given time. It was 
shown that the condition under which this is a good approximation is that 
the square of the ratio of the impact velocity to the speed of sound in 
the gas he small compared to unity. Fortuitously, the limitation on impact 
velocity due to skin strength of practical materials automatically satis- 
fies this condition if inflating gases with high speeds of sound, such as 
helium or hydrogen, are used. 

The preliminary analysis of reference 1 has left some questions 
unanswered. In reference 1 the mass of the suspension cords was assumed 
either to he negligible or to he taken into account approximately hy being 
included in the mass of the package at the center. The present report is 
a design study of the inflated sphere landing vehicle in which definite 
consideration is given to the mass of cords required. The required mass of 
cords depends, of course, on their mass density and strength; and this mass 
is found to he of very significant magnitude for constant-area cords 
constructed of practical materials . 1 There is also the question of how 
critical will he the effects of deviations from the design impact velocity 
and design initial gas temperature, and the effects of errors in the time 
at which the skin is ruptured. A knowledge of these effects would aid in 
the selection of the design parameters. A principal function of this 
report, therefore, is to consider these effects and present the results. 

The results of the analysis show essentially the behavior of the vehicle 
when deviations from design conditions occur at impact; hence indications 
are given as to how the design parameters should he chosen for reliability. 
The actual tolerances in the deviations from design conditions are, in fact, 
found. 


1 Further investigation has shown that if the suspension cords were 
tapered so that a constant stress were present over the length of one 
cord at a given time, the payload mass which could he landed, for some 
design conditions, would he significantly higher than for the case of 
constant-area cords. (See appendix.) 
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SYMBOLS 

acceleration of the sphere center 

a, dimensionless acceleration 
cross-sectional area of one suspension cord 
specific heat at constant volume 

effective longitudinal modulus of elasticity of each cord 
cp 

— , distribution of force in the cords 
^o 

resultant vertical force due to the gas pressure 

gravitational acceleration on the earth 

drop altitude above the moon surface 

total mass of the landing system, including payload 

mass of a particular part of the system, depending on the 
subscript 

maximum acceleration in earth g T s 

total number of suspension cords 

pressure of the inflating gas 

pressure of the atmosphere outside the vehicle 
Pa 

Pi 

radius of the sphere 
gas constant 

time measured from first instant of impact 



temperature of the inflating gas 

velocity of the sphere center, ^ 
u 

I u x | 
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distance from the impact surface to the sphere center 
r 

arc cos (-y) (see fig. l) 
given by equations ( 28 ) and ( 36 ) 
given by equations (4l) through (4-3) 
ratio of specific heats 

vertical distance the package deflects as a result of the elastic 
extension of all cords 

fractional error in initial temperature 

fractional error in the initial impact velocity 

fractional error in the initial energy ratio 

critical value of e*; the value for which y 2A = 0 

radial coordinate (see fig. l) 

angle indicated in figure 1 

defined by equation ( 20 ); may vary from l/3 to l /2 

mgCyTi 

initial energy ratio, — 

mass density 
stress 

tensile force at t| = r in each cord at angle 0 
total vertical force on the skin from the cords 



Subscripts 

0 value at 0=0 

1 condition at first instant of impact 

2 condition at maximum compression, when the velocity is zero 
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A value of a parameter in the actual case 

c suspension cords 

D design value of a parameter 

f condition at final impact; that is either at y = 0 or at end 

of rebound 

g inflat ing ga s 

p payload package 

r condition at the time skin is ruptured 

s skin 

v portion in motion with velocity u 

ANALYSIS 

Properties of the Sphere During Impact 


As in the analysis in reference 1 ; it is assumed that the skin material 
is flexible, but nonstretchable, and thus that the part of the sphere skin 
not in contact with the impact surface retains its spherical shape, because 
it is supported by the internal gas pressure (see fig* l) . Then each 
particle of that part of the skin moves with velocity u, which is also 
the velocity of the center of the sphere. If it is also assumed that the 
suspension cords are flexible but nonstretchable, then the suspended pay- 
load package remains coincident with the sphere center and moves with 
velocity u, and the parts of the cords not in contact with the stationary 
part of the skin which is held on the surface may also be assumed to be 

moving with velocity u. Thus, all particles of the system which have 

struck the impact surface are assumed to have velocity zero, and all others, 

except for gas particles, are assumed to have velocity u. At the first 

instant of impact, t = t x , all particles of the system have velocity u^ 

At the instant of maximum compression, t = t 2 , all particles have velocity 
zero. Also at t = t 2 , if the sphere has compressed nearly to a hemisphere, 
the payload has just reached the surface at the instant its velocity is 
zero, and the maximum permissible amount of energy has been absorbed by 
the inflating gas. 
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Equation of Motion 


The total mass of the system, m, is considered to consist of four 
parts: 


m = mp+m s +m c +mg (l) 

where irip is the mass of the payload package suspended in the center, mg 
is the mass of the skin, is the mass of the cords, and mg is the 

mass of the inflating gas. One may follow a procedure similar to that 
used in reference 1 and define that portion of the mass which remains in 
motion, by 


% — m SV + ®CV + (2) 

where m sv and mcv are, respectively, the parts of the mass of skin and 
mass of cords which are in motion with velocity u. 


It was shown in reference 1 that, if the gravitational force and 
flexural stresses at the impact surface are neglected and if the uniform 
gas approximation is used for the internal gas pressure, application of 
Newton’s second law of motion to the mass m^ leads to the differential 
equation of motion: 


1 d(u 2 ) 3(7 - 1)1 m _2 ,[/ 4 Y 

a = 2 — = 8 5* (1 ' y } LU + 3 y - r) 



(3) 


where p a is the outside aerodynamic pressure (taken to he 
constant at the value of the local atmospheric pressure), £ 
energy ratio defined as in reference 1 hy 

mgCyT! 

1 (l/2)mu! 2 

and y = r y is the distance from the impact surface to the 
Equation (3) applies to the present analysis if the relation for my 
given hy equation (2) replaces that used in reference 1. 

The differential equation, (3), can he solved if my/m is written as 
a function of y. As shown in reference 1, m sv is given hy 

m sv = - m s (l + y) (5) 

The portion of the mass of cords which has not struck the surface, m cv , 
may he found as follows: let tj he the radial coordinate and G he the 

cone angle as indicated in figure 1. The mass density and cross-sectional 
area of one cord are p c and A c , respectively; hoth are constants. The 


essentially 
is the initial 


(4) 


sphere center . 
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mass of an element of length of one cord, dr), is then p c A c dT). If the 
total number of suspension -cords is N c , then the number of cords in the 
angle d 9 is 

dN c = ^ N c sin 9 d 9 (6) 


Equation (6) is, of course, approximate if N c is finite. Therefore 
the mass in a radial element, dr], and an element of 0, d0, is 


1 

dm c = “ N c p c A c sin 9 d 9 dq 


(7) 


but the total mass of the cords, if A c is constant, is 


me — N c p c A c r 


Equation (7) ma y thus be written 

dm c = — sin 9 d 9 dT) 

2r 

Then m cv is 

p9=a nr\=r p9=x pj)~y/cos 9 

m cv = / &&C+ / toe 

J 0=0 J T)=0 J 9=a J T]=0 


(8) 

(9) 


(io) 


where 


a = arc cos(-y) 


(id 


Substitution of equation (9) into (10) and use of the result, along with 
equation (5) > in equation (2) give the following relation for m v /m as 
a function of y: 


nfy nip m s m c + y - 

"m" - m” "m" \ 2 / m \ 2 

For specified values for the masses of components of the system, the 
differential equation, (3), may now be integrated to obtain u as a 
function of y. Simultaneously, the equation 

dt/dy = l/u (13) 

may be integrated to obtain y as a function of t . The initial 
conditions for equations (3) and (13) are: 

t = 0, y = l> u = -1 (l^-) 


y in y 


( 12 ) 
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The value of £ to he used in equation (3) may either he specified 
or he found, for a desired value of y 2 (the value of y when u = 0 and 
t = t 2 , the instant of maximum compression of the gas), hy the procedure 
discussed in reference 1. Thus the motion of the sphere during impact 
is determined if mp/m, m c /m, m s /m, and either £ (defined hy eq. (4)) or 
y 2 are specified. 


Performance Requirements and Capabilities 


In addition to determining the motion, that is, acceleration, velocity, 
and location histories of the sphere center, equation ( 3 ) can also he used 
in conjunction with certain other relationships among the various param- 
eters to determine the performance requirements and capabilities of the 
landing system. The performance includes mass of skin required, mass of 
gas required, mass of suspension cords required, mass of payload package 
allowed, and the required radius of the sphere. It was seen in reference 
1 that these characteristics depend on initial gas temperature, gas 
constant and ratio of specific heats in the gas, a stress to mass ratio 
in the skin, impact velocity, final volume ratio, and outside atmospheric 
pressure. The radius required was also seen to he inversely proportional 
to the allowable maximum acceleration of the payload. Because the mass 
of suspension cords is included in the present analysis, the performance 
characteristics will also he a function of the strength to mass ratio of 
the cords . 


The procedure to he used in evaluating the performance is as follows: 
Specify values for the parameters C7 sl /p s , cr c 2 /p c , 7 , R, T x , p a , y 2 , and 
mp/m. Then develop four equations for determining the proper values of 
m c /m, m s /m, mg/m, and nr in relation to each other and to the specified 
parameters. Integrate the differential equation, (3)> to determine £, 
from which can he found the corresponding impact velocity ui- 


One of the four required relationships is the same as equation ( 33 ) 
of reference 1 , an expression relating the mass of skin to the mass of gas 


= 3RTi(l-p a ) 
m 2 (a sl /p s ) m 


(15) 


where cr s i is the stress in the skin before impact . 
from equation (l) , is 

mg mp m s m c 

m m m m 


Another relation, 


( 16 ) 


The mass of suspension cords required can he found in terms of the 
maximum stress in the top cord (which can he chosen to he less than the 
ultimate strength of the cord material) . 
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The cords can support only radial tensile forces . Let the distribution 
of force in the cords he given by 

9 = cpof (17) 

where cp 0 , the force in the vertical cord (9 = 0) , is a function of y and 
f is a function of y and 9. The vertical force from the cords on the 
skin in the angle d0 is 


d4> = 9 cos 9 dN c 

= ^ N c 9 0 f cos 9 sin 9 6.9 (18) 

where use has been made of equations (6) and (17) • Only the cords from 
9 = 0 to 9 = rt/2 are in tension during the impact. Therefore, the total 
vertical force transmitted by the cords to the skin is 


0 = 


9=n/2 

9=o 


d<£ = g NcTbd 


(19) 


where 


vjt / 2 


d = 


f sin 9 cos 9 69 


( 20 ) 


The total vertical force exerted by the gas pressure in decelerating 
mv is, from equation (8) in reference 1, 

Fp = m v (du/dt) (21) 

Because it is assumed that all parts of m v , that is, m sv , 1%^, and nip, 
have the same velocity, u, the rate of change of momentum to m sv is 
m sv (du/dt) . The remainder of F p is transmitted to the mass (mp + n^y) . 
Therefore, the rate of change of momentum transmitted to (mp + 1%^) is 
(mp + m cV ) (du/dt) and is equal to the force transmitted through the cords: 

<5= (mp +mcv)(du/dt) (22) 

Thus, from equations (17)^ (l 9) > and (22), the force at q = r in one 
cord at angle 9 is 

^ = (m P + mcv) at f (23) 

The maximum force at a given time during the impact, as will be seen, is 
in the vertical cord at 9 = 0, t] = r, and the maximum occurs at the 
instant of maximum compression (t = t 2 ) . (When y 2 = 0, the maximum 
occurs a very short time before t 2 but is always nearly identical to the 
value at t = t 2 .) Therefore, 




r 
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Vx 

The maximum acceleration is 


2 

N c H 2 


(mp + mcv 2 )(% 


(2k) 


(du/dt) 2 = ng e (25) 

that is, n is the number of earth g f s maximum acceleration. The left 
side of equation (24) may he written as 

q^nax = a c2A c (26) 

where cr c2 is the maximum stress in the vertical cord. Substitution of 
equations (25) and (26) along with equation (8) into equation (24) then 
gives 


te(cr c ^/pc) 
m ~ m L 2g e (nr) 



(27) 


where, from equations (9) and (10), 

P 2 = = \ (1 + 72 - y^Zn y 2 ) ( 28 ) 

Equation (27) supplies the third of the four required equations relating 
the various parameters . 


It is now seen that the required m^/m depends on the strength to 
mass ratio of the cords, o^/pc* 14 also depends on the value of \± 2 , 
given by equation (20) when f = f(y2> 0) . It is noted that exact knowl- 
edge of f, the distribution of force in the cords, is not of critical 
importance in this study because only [i 2 , the integral of the product of 
f times a function of G over G from zero to jt/2 at y = y 2 , is 
involved in the final results. The final results can be given in terms 
of [1^2(0*02 / p c )]* Then \± 2 can be estimated, and the effect of the choice 
of [± 2 can be considered to be equivalent to a change in the allowable 
value of o^/pc* Limits to the variation of f(y 2 , 0) > and thus to the 
choice of \± 2 , can be shown. In those cases for which the tensile stress 
at all points in the upper hemisphere remains greater than zero, the 
shape of the skin will necessarily remain exactly spherical, regardless 
of the forces acting normal to the skin, because of the assumption of high 
modulus of elasticity in the skin. The skin above a horizontal diameter 
will then act as a rigid shell as long as cr s > 0 at all points in the 
upper hemisphere. Then, if the mass of one cord is negligible in compar- 
ison to mp and if (Sp/r) « 1, cp is given by 

AL = &pCos 9 = cpL/AcE , (0 <0< it/2) (29) 

(where AL is the small change in length of each cord whose original 
length is L = r, the sphere radius), as long as ct s > 0 at al 1 points in 
(0 < 9 < it/2) . Thus if cr s > 0 
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f = cpApQ = cos 0 , (0 < 0 < at/2) (30) 

and, from equation (20), if cr s > 0 at all points in 0 < 0 < x /2, 

n = l/3 (31) 

For all cases the skin stress is everywhere greater than zero in the 
upper hemisphere at least for a time at the beginning of the impact. In 
some cases this condition will exist until t = t 2 , and thus p 2 will 
have the value 1 / 3 . For other cases , because the maximum force per unit 

area on the skin is at 0 = 0 (eq. ( 30 )), the skin stress at 0 = 0 may 

be reduced to zero at some time during the impact (as the acceleration 
increases, and thus as the forces in the cords increase). After that 
time, the cosine cord-force distribution given by equation ( 30 ) is no 
longer valid because the upper hemisphere can no longer be considered 
entirely a rigid shell. If f were required to maintain the cosine dis- 
tribution as the acceleration increased further, (cr s ) 0 would tend to 
become negative and thus a small part of the skin at 0=0 would tend to 
collapse because it is flexible. However, if this happened, the force in 
the vertical cord would tend toward zero while the cords at higher 0 
carried the increased load. But, if the force in any of the cords became 
reduced because of skin collapse, the cause of the reduction in the skin 
stress would be removed. Therefore the skin should tend not to collapse, 
and the cord force distribution would simply accommodate accordingly; 
that is, f becomes nearly constant in the region of zero skin stress. 

It is seen that, for those cases in which mp/m and iric/m are sufficiently 
large, large cord forces are present during the impact, and a region of 
zero stress may occur; hence the function f becomes implicitly involved 
with the stresses in the sphere skin. A detailed analysis of the skin 
stress distribution would be required to determine f(y 2 , 0), and thus 
P 2 , exactly. 


As noted above, for purposes of obtaining results in the present study, 
it is sufficient to determine limits to the estimated values of p 2 - It 
is therefore noted that, after (cr s ) 0 becomes zero, the cord force dis- 
tribution accommodates so that the cord force is nearly constant in the 
region of zero skin stress. The limiting case is therefore a constant 

value of cp over 0, that is, f = 1, for which p = l/2. Thus p 2 will 

be as low as 1/3 in some cases, but, for those cases where mp/m and m c /m 

are high, it may be nearly l/2, and the effect of varying p 2 is simply 

a change in the value of the allowable cr c2 / p c . 


The final expression of the four required (as discussed above) may be 
obtained from equation (3) "by noting that the maximum acceleration is 
given by 

%iax = a 2 = (r/u 1 2 )ng e (32) 


and by making use also of equation (4) and the equation of state of the 
gas at t = t x : 



mgRT x = mg(/ - l)c v Tx 


(33) 
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The result, which can he used as the fourth required equation relating 
the parameters, is 


3 RTi ^ m / — 2 \ 

nr = 1 ~ C 1 “ 72 ) 

4 g e m 


LV2 + 3y 2 - y 2 3 


“ P« 


(3M 


where, from equations (12) and (28): 


m v2 _ , m s + y2^\ , n me 

m m m \ 2 J + ~ 


( 35 ) 


Equation ( 34 ) may now he combined with equations (15), (l6), and (27) 
to obtain independent explicit relations for m s /m, mg/m, nic/m, and nr 
in terms of mp/m and the various other parameters involved. For con- 
venience in the ensuing manipulations , the following parameters are 
defined: 


Pi = M2 ( h C2 / p c ) / 2g e 

P2 given hy equation (28) 

P 3 = ( 3 / 2 )RTa. 

34 = ffsi/ p s 

P = = Hs4(r + 37 ,- yj - 1 

3 e = \ (l + Is) 


(36a) 

(36b) 

(36c) 

( 36 d) 

(36e) 


Substitution of equations (l6) and ( 35 ) into ( 34 ) gives 


nr 


P 3 3 5 (l - mp - m s - me) 
Sp + 3 6 m s + P2™c 


(37) 


where m( ) indicates m( )/m. Equations ( 15 ) and (l6) can be combined 
to eliminaue mg/m and to obtain 


P 3 (l - p a )(l - n*p - me) 
P 4 + P 3 (l - Pa) 


( 38 ) 


Substitution of equation (38) into ( 37 ) then gives an expression for the 
product nr in which the only mass ratios involved are nip and rr^: 
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P3P4Ps(l - mp - me) 

[(B 4 + P 3 (l - p a )](mp + P 2 m c ) + P 3 Ps(l - Pa)(l - % - me) 


(39) 

Equation (39) may then he substituted back into equation ( 27 ), from which 
a quadratic expression for tUq is obtained. The result is 



where 

P 7 = PiP2tP4 + ^ 3(1 ~ Pa) ] — PiP3Pe(l - p a ) + P 2 P 3 P 4 P 5 (^-l) 

38 = Pit 34 + 33(1 - P a )]®p + 3l333s(l - P a ) (l - mp) 

- 32333435(1 - mp) + 333435mp (^ 2 ^ 

P 9 = 3 3 p43 5 mp(l - mp) (^3) 

The positive sign on the square root must be used to give positive values 
of mc/m. As Hp approaches zero, me approaches a positive finite 
value. Thus, if there are cords in the vehicle and if nL, = 0, the plus 
sign must be used, and a certain mass of cords is required. However, if 

there are no cords in the vehicle, = 0 and the minus sign may be used 

to give me = 0 (since P 8 is negative). 

If mp/m is given, the corresponding required mc/m can be found 
directly from equation (4o) . Then the required value of m s /m is found 
from equation ( 38 ), and the required mg/m, from equation (l 6 ) . Finally, 
the product nr is found from equation (37) • Each of these results, 
however, corresponds to a certain impact velocity, u x . The impact 
velocity can only be found by integrating the differential equation, (3), 
and using the procedure outlined in reference 1 to obtain the correct 
value of | corresponding to the desired y 2 and the other input condi- 
tions . The correct value of u x is then found from equation (4), which 
can be rearranged in the form 


ui = - 


^3 3 mg 


3(7 - 1)1 


(44) 


The drop altitude above the moon surface corresponding to a given impact 
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velocity is 


hm = Ux 2 / 2g m = u x 2 / (10.94 ft/sec 2 ) ( 45 ) 

where g m is the gravitational acceleration on the moon. For convenience 
this correspondence is shown in figure 2. 


Effects of Deviations in Gas Temperature and Impact 
Velocity From Design Conditions 


In the design of the inflated sphere landing vehicle a certain initial 
gas temperature, T x , and a certain impact velocity, uj., would he assumed. 
The appropriate design mass ratios and sphere radius would then he calcu- 
lated from the equations in the preceding section. In the actual landing, 
however, the gas temperature and intact velocity may differ slightly from 
the design conditions . The effects of these deviations may he inconse- 
quential or very critical, depending on whether a certain condition was 
overestimated or underestimated. The choice of the design parameters will 
thus he influenced hy the effects of deviations from design conditions. 

For example, suppose the landing vehicle is designed for a given u x 
and Tx so that it compresses to a hemisphere, that is, y 2D = 0. Then, if 
the actual impact velocity is slightly higher than the design u x , the 
payload will strike the surface with some residual velocity because it has 
not been completely decelerated to zero velocity at y 2 = 0. However, if 
the actual velocity is slightly less than the design Ux, then the maximum 
compression will occur at y 2 slightly greater than zero; that is, the 
payload package will come to rest a short distance above the surface, then 
drop to the surface when the skin is ruptured. The effect of the error 
in ux in this case would he much more tolerable than in the former case . 
The design y 2 should therefore he chosen so that there is a high proba- 
bility that the latter case will occur rather than the former. 

Errors in the initial gas temperature will have similar effects. A 
negative error (value below design) in Tx will cause y 2A to he less 
than y 2D , or the payload will strike the surface with some residual 
velocity, whereas a positive error in Tx will have the effect of decel- 
erating the payload to zero velocity at y 2A greater than y 2 p because 
of the higher gas pressure. In the case of positive error in Tx, the 
effect may not necessarily he inconsequential because the final accelera- 
tion will he higher than the design value. 

It is the purpose here to estimate the magnitudes of the effects of 
deviations from design conditions and hence to show how the design 
parameters should he chosen for reliability. 

The landing vehicle is designed for the conditions T-^ and Uj^, as 
well as the desired y 2 j)> n^, and other parameters involved. The 
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corresponding design values of m s /m, mg/m, nic/m, mp/m, r, and {■ are 
found "by application of the theory in the previous section so that u = 0 
when y = where ^ = u /| u 1dI* In the actual case, in order for the 

initial condition, equation (l4), to apply in the integration of equation 
(3), u must he considered to he u/|u 1A |, hut this difference does not 

affect the integrated results of equation (3). The only parameter in 
equation (3) which can differ from the design condition is the initial 
energy ratio, g; that is, | A may not he equal to and hence y 2A 

will not equal y 2 p. An error in | may occur because of an error in 
T]., Ua, or hoth. 

Denote a fractional error in g hy €|j that is, 

t A = 0- + egU D C 46 ) 

Also denote errors in T x and u x by and € u , respectively. Thus 

Tp A = (1 + € T )t 1d ( 4 t) 

u ia - (l + e u^ u lD 

Therefore, since 

g = 2(m g /m)c v (T 1 /u 1 2 ) ( 4 9) 

the result from equations (46), (47), (48) and (49) is 


1 + e 


I “ 


1 + Giji 

(1 + € U ) 2 


(50) 


For no error in velocity^ then^ the error in | due to an error in 
temperature is simply 


( e |)e u =0 " e T 


(51) 


and for the case of no error in temperature, the error in g due to an 
error in velocity is 

. _ ~ e u( 2 + e u) 

(s |h^o - (1 + eu)2 

In the cases where e u is small , 



( £ |>€ T -0 
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For the cases where is positive (e.g., e T > 0 and e u =0, or e u < 0 

and e T = 0) the actual 5 value of y 2 is greater than the design value, 
that is, for > 0, 

y 2A > y 2D m 

For the cases in which is negative (e.g., e T < 0 and e u = 0, or 
e u > 0 and e T = 0) either the payload strikes the surface with final 
residual velocity^ uf or the payload does not reach the surface, but the 
actual value of y 2 is less than the design value; that is, for e» < 0: 
either * 


or 


u 


= u^ when y = 

0 - y 2A < y 2 D 


0 


(55a) 

(55b) 


The effects of deviations from the design conditions for a given case 
(given T 1d , y 2D , mp/m, and other input parameters and the calculated 
values of !p, uqp, m s /m, mg/m, and m c /m) are thus found by integrating 
equation (3) again, using the same values for the mass ratios, but using 
in place of £-p. The results of the integration (y vs. a, u, and T) 
dhow either the value of y 2A w ben u = 0 (and the corresponding value of 
a PA ) or the value of up when y = 0. It should be noted that 
Uf = Uf/ | u 1A | and therefore that 

Uf = (l + e u )|u 1D |uf (56) 


The ratio of actual maximum acceleration to design maximum acceleration, 
for those cases in which the payload does not strike the surface, can be 
found as follows: 


From equation (25), 


Therefore, 


(u x 2 /r)a 2 = ng e 

(57) 

a 2 D = g e n D r / u lD 

(58) 

a 2A = g e n A r/u 1A 2 

(59) 

n A / n D = ( a 2A / a 2 D^-*- + e u) 2 

(60) 


and thus. 


Effect of Error In the Time at Which Skin Is Ruptured 


It is desired to have the skin rupture at time t2p i when the actual 
velocity is zero. It is assumed that some method is used to rupture the 
skin instantaneously around a great circle of the sphere . Each half of 
the sphere skin 'will remain attached hy half the suspension cords to the 
payload. The gas will escape simultaneously in all directions outward 
from the center in the plane of the rupture. The vehicle will tend, 
therefore, not to he propelled hy the escaping gas. 

Let t r denote the actual time of skin rupture, which may not he the 
same as the most desirable time of skin rupture, t2A* ^ At is the 
time error, then 

t r = t2A + At (6l) 

where At may he positive or negative and where, it must he noted. 


At = — - — At = 

I U 1A | 

to satisfy condition (lU) on equations ( 3 ) and ( 13 ) • 

If the error At is negative, the inflating gas will be released 
from the vehicle before the payload has been decelerated to zero velocity. 
The decelerating pressure force will therefore be eliminated, and, if 
this occurs instantaneously, the payload package will strike the surface 
with the velocity it had at the time of skin rupture; that is, u^> = u r . 

The velocity u r can be found directly from the integrated results of 
y versus a, u, and T, provided the velocity u-^ is used in the conversion 
from dimensionless to dimensional parameters. 

If the error At is positive, the sphere will reach its maximum 
compression and begin to rebound before the skin is ruptured. The magni- 
tude of the velocity will be essentially the same at t = ^2k + ^ as 
at time t 2A - At, but the motion is in the opposite direction. This 
velocity, u r , however, is not the same as the final impact velocity of the 
payload after it falls back to the surface. The reason for this is that, 
at t = t r , nearly half of the mass of the skin and cords is at rest; 
therefore, some of the rebound momentum is taken up in putting this 
stationary material back into motion. The velocity of the payload at 
final impact, uf, is therefore nearly the same (in magnitude) as the 
velocity after the entire system has just left the surface. It is assumed 
here, of course, that u r is great enough that the sphere rebounds com- 
pletely and with a final velocity great enough that the increment in 
velocity on the final impact due to the difference in height (equal to the 
sphere radius) is negligible. Equating the momentum at t = t r to the 
momentum at the instant the entire mass has left the surface during the 


(l + e„) luin 1 


At 


(6 2 ) 
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rebound, assuming the gas is moving at the same velocity as the rest of 
the system prior to t r and is dispersed simultaneously in all directions 
at t = t r , one finds 

(nVp + mg)u r = (nip + m s + m c )(-u f ) 


Thus 


( u r 

u f = - u ia(— 


1A 


(m V r/ m ) + (mg/m) 
-(mp/m) + (m s /m) + (mc/m)- 


(63) 


where mvr/m is the value of i%/m given by equation (12) when y has 
the value corresponding to t = t r . Of course, t r is given by 



(610 


where At is positive in this case. 


DISCUSSION OF RESULTS 


Performance Requirements and Capabilities 


In the analysis of the performance of the inflated sphere landing 
vehicle, the mass of the suspension cords has been considered. The results 
are shown in figures 3 through 7 . The design parameters, including the 
distribution of the component masses of the system and the radius required, 
are shown plotted versus impact velocity (the corresponding drop altitude 
above the moon is shown by fig. 2). The results are given for the case 
of no outside atmospheric pressure and for hydrogen as the inflating gas 
with 7 = l.ll-l and R = 766.5 ft-lb/lb m °R. The skin strength parameter 
for all results shown is cr sl /p s = 10 6 ft 2 / sec 2 , a reasonably attainable 
value, as is discussed in reference 1. Several values of ypp and 
various combinations of the values of initial gas temperature of 200° and 
500° R and the values of the product of p 2 and the strength-mass ratio 
of the cords of 1.2x10 s and 1.6x10 s ft 2 /see 2 are employed. These values 
were computed assuming p 2 = 0A and using for (cr c2 /p c ) typical values 
of yield strength to mass density ratio of either nylon cords or steel 
wires: 3X10 6 ft 2 /sec 2 and ^X10 6 ft 2 /sec 2 . Therefore, if actually 

p 2 = l/2, then the results computed for [Ms(cTc2/pc) ] = 1-6X10 6 ft 2 / sec 2 
would correspond to Gcz/ Pc = 3-2X1O 0 ft 2 / sec 2 , or if \i 2 actually is 
l/3 j then the results computed for [fe(cr c2 /p c ) ] = 1* 6x10 s ft 2 / sec 2 
would correspond to cr c2 /p c = 8x10 s ft 2 / sec 2 . The value of p may 

remain at l/3 during the impact in some cases, whereas in others it may 
vary to a number between 1/3 and l/2. In those cases where mp/m and 
mo/m are large, p 2 niay be nearly l/2. It should be noted that, in the 
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calculations for figures 3 through 7, there has been no consideration of 
mass of any devices for attaching the cords to the sphere skin. It is 
presumed that a method of attachment could be developed which would not 
require a significant extra mass at the cord ends, but which would never- 
theless minimize stress concentrations in the sphere skin at points of 
attachment by distributing the force from a cord over a circle or an area 
rather than applying it at a single point. 

The performance results given by figures 3 to 7 can be used as follows: 
For a given desired impact velocity, the ratio of the mass of the payload 
package to the total mass of the system is found from figure 3 - (Note 
that these results depend strongly on the impact velocity. Their depend- 
ence on the allowable strength-mass ratio of the cords is also significant, 
just as it was found in reference 1 that the results depend strongly on 
the skin strength parameter, cr s i/Ps*) ^ required mass of skin, mass of 
suspension cords, mass of gas, and value of the product nr are found from 
figures 4 through J. For a given allowable acceleration, then, the radius 
required is known. For the desired total mass, m, the required skin 
thickness may be calculated from 


m s = 4Ttr 2 p s 8 

using the value of m s /m from figure 4 . The quantity N c p c A c can then 
be found from equation (8) using the value of nic/m obtained from 
figure 5 . Therefore, if a given cross-sectional area for the cords and 
the cord material are chosen, the number of cords required is determined. 
(As a check to insure that the volume of cords is not a significant per- 
centage of the sphere volume so as to affect the results, one can calculate 
the ratio 


mc/Pc 
( 4 / 3 ) hr 3 

For practical cases, this ratio is very small. In a numerical example 
given later in the discussion, if the total mass is less than 1300 pounds 
and if steel wires are used as the cords, this ratio is less than 10 -3 .) 
Finally, the ratio of mass of gas to total mass is found from figure 6 . 
The internal gas pressure may then be determined by the equation of 
state, (33)* 

Figure 8 has been plotted to illustrate the comparison with respect 
to temperature for the ratio of the payload package mass to the total 
mass of the system. It is seen that mp/m varies little with varying 
initial tenqperature of the inflating gas. 

In figure 9 is given a comparison of the present results with those 
of reference 1 . It is seen that the results for mp/m differ signifi- 
cantly, because the required mass of cords has been considered separately 


= 1 / m 
4 jt ' p c r 3 


me 

m 
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from the payload mass in the present report, but that they very nearly 
agree if the value of mp/m in reference 1 is assumed to include both 
payload and required cord mass. (The reason that (mp + m^/m / (mp/m) re f . q 
is that the center of mass of me V does not remain coincident with 
that of mp.) Note that at the higher velocities the mass of cords 
required in the present analysis would exceed the value of nu given by 
the analysis of reference 1. Also note that the mass of cords required 
(for constant-area cords) does not approach zero as mp/m approaches 
zero. (See also fig. 5.) 


Effects of Deviations in Impact Velocity 
and Initial Temperature 


As indicated in the Analysis, e g is a convenient parameter for 
relating the deviations in initial temperature and impact velocity to the 
resulting effects of such deviations. The effects are therefore shown 
in subsequent figures only as effects corresponding to deviations in g. 

For convenience, figure 10 shows a graphic representation of the relation- 
ship between the deviation in | and the corresponding deviations in ui 
and Ti plotted directly from equation (50 ) . It can be observed from 
these curves that the largest errors in £ occur for positive errors in 
Tx and negative errors in Uq. However, this condition gives a positive 
error in g, and, as mentioned 'in the Analysis, this is not the critical 
case . 

The effects of deviations in uq and Tq, through the relationship of 
e u an & £rp , with €| given in figure 10, are presented in figures 11, 12, 
and 13- In those cases where eg is negative and where the payload has 
not been decelerated to zero velocity at y = 0, the residual velocity, 

is given by figure 11. The values of uqp correspond with the design 
parameters given by figures 3 through 7- Curves are shown for various 
values of eg and various y 2 D' The results are left in terms of e u 
so that only eg need be specified. In part (a) of figure 11 the results 
are for Tqp = 500° R, whereas those of part (b) are for T 1D = 200° R. 

One may note that the results shown are very nearly identical for this 
temperature range except for the upper end of the curve for eg = -0.40 and 
y 2 p =0. It is seen that, when negative deviations in g occur, the pay- 
load may strike the surface with a high velocity if the design value of 
y 2 is too low. As is increased, this is much less likely to occur. 

In fact, for y^^ =0.2, the payload does not strike the surface even for 
eg = -0.40 under the design conditions indicated in figures 11(a) and 11(b). 

In the cases where the payload does not strike the surface before the 
vehicle rebounds, the velocity will be zero at y = y c A and the maximum 
acceleration will be n^ earth g’s. This occurs whenever eg is positive 
and in some cases when eg is negative. Curves of y 2 A versus ” U 1D 
are given in figure 12. 
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A result of greater interest than the value of ^ or ^ ese cases 

is the ratio of actual to design maximum, acceleration, shown in figure 13 • 
Each point in figure 13 corresponds to a point in figure 12. It is noted 
that, for many cases, the actual maximum acceleration exceeds the design 
value. However, n^/n^ is, in general, not much greater than unity. This 
ratio tends to he near or less than unity especially at the higher veloci- 
ties, which are of greater practical interest. The higher values of 
(n^/n^)/(l + e u ) 2 occur for positive values of eg and for low design 
impact velocities. It is seen from figure 10 that large positive values 
of eg occur only for negative values of e u . Thus the actual value of 
n A/np will he less than (n A /n;p)/(l + e u ) 2 plotted on figure 13 when eg 
has a large positive value. Figures 13(a) through 13(d) show that 
is slightly greater for higher J2 D when is positive, hut that 

n /\/ n D l°' w ’ er i* or higher y2D when eg is negative. The results for 
n l/ n D are use i’ u i n °t only i n determining the actual maximum acceleration 
of the payload hut also in finding the actual maximum stress in the cords. 
For this purpose, it can he shown from equations (24), (25), and (26) that 

°c2A _ V 

a c2D n D 

where (3 2 A is the function of y 2 A given by equation (28) and P 2D 
is the same function of y 2 p # 

The critical deviation in |, eg*, may he defined as the fractional 
error, for a given design case, which causes the sphere to just compress 
to a hemisphere; that is, the deviation in £ has a value such that 
y 2A =0. If y2D = ^en, of course, €|* = 0. Plotted in figure l4 are 
the values of -€£* corresponding to design impact velocity. It will 
again he noted that all cases considered here are calculated for a vehicle 
designed according to figures 3 through 7- It is seen that if the vehicle 
is designed with y 2 p only slightly greater than zero the critical value 
of £ g is fairly large. For example, a vehicle designed for = 0*^0 > 

an impact velocity of 500 ft/sec, and the other conditions as indicated 
in figure l4(a), can experience a negative deviation in £ as large as 
25 percent and not have the payload strike the surface. If, instead, 

0.20, then the negative deviation in £ can he as large as 45 
percent . It may he noted that these results depend very little on T^p 
in the range shown here . 

As an example of the use of the above results, suppose a moon-landing 
vehicle is designed for an impact velocity of 500 ft/sec (which corresponds 
to a drop from an altitude of 23,000 ft above the moon f s surface), an 
initial hydrogen gas temperature of 500° R, an allowable skin strength 
parameter cr sl /ib s = 10 6 ft 2 /sec 2 , allowable strength-mass ratio of the 
cords of 4xl0 6 ft 2 /sec 2 (assuming \± 2 = 0.4), and y^p = 0.10. Suppose 
also that the maximum acceleration allowed is 2000 earth g T s. From the (a) 
parts of figures 3 through 7, the results are: mp/m =• 0.44l, m s /m = 0 . 258 , 

ir^/m = 0.287, mg/m = 0.0l4, and nr = 11,200 feet. The radius required is 


(mp/m c ) + 

+ Bo-n- 


(65) 
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therefore 5-6 feet. Because the vehicle is designed for y^D = 0.10, 
there can he some error in £ and the payload will still not strike the 
surface (if the size of the payload itself is not considered). The 
critical error, from figure l4(a), is €g* = -0.25. This might correspond, 
for example, from figure 10 (h), to e u = 0.20 and e T = 0 . 08 ; that is , the 
actual impact velocity can he 20 percent higher than design, or 600 ft/sec, 
if the temperature is 8 percent higher than design, or 540° R, before 
the payload will strike the surface with any residual final velocity. 
Suppose, however, that the deviation in £ exceeds the critical value 
(negatively) . For example, suppose that the initial gas temperature is 
13 percent lower than the design temperature, or only ^35° R, and that the 
impact velocity is 20 percent high; then e * = -0.40. One can see then, 
from figure ll(a), that -Uf/(l + e u ) = 222 ft/sec and therefore that the 
payload will strike the surface with a final velocity of 2 66 ft/sec. 

Consider now the opposite case in which the deviation in | from the 
design condition is less than the critical value. Take, for example, 
c ^ = -0.10. This deviation in £ would occur if there were no error in 
impact velocity, hut a -10 percent deviation in temperature, that is, 

Tia = 450° R. It is found from figures 12(h) and 13(h) that the payload 
then reaches zero velocity at y>>A = 0.064 and that the maximum 
acceleration, in earth g f s, is n A = 0 . 992 ( 2000 ) = 1980 . 

Considering the above results, one can now evolve a "design philosophy" 
for the inflated sphere landing vehicle; that is, one can determine how 
the design parameters should he chosen for reliability. From figure 3 it 
is seen that the payload is somewhat decreased when J 2 J) i s increased 
above zero . For small increases in the decrease in payload is small . 

Figures 7(a) through (d) show that the radius required is not significantly 
affected by increasing y 2 D # Figures l4(a) through (d) then show that the 
critical (allowable) deviation in £ is greatly increased (in magnitude) 
by increasing y 2 p slightly. If the critical deviation in £ is only 
slightly exceeded negatively, figure 11 shows that the final impact 
velocity of the payload is very high, and therefore critical. Figure 13 
shows that the consequences of negative deviations in £ less in magnitude 
than the critical value and of positive deviations in £ are insignificant. 

It is now apparent that a value of y2D shou ld used which is as 
high as possible but which allows a sufficient payload-landing capability. 
The actual y^p used will, of course, be dictated to some extent by the 
expected deviations from design conditions. It should be noted that the 
choice of design y 2 should also take into account allowances for the 
size of the payload, for slight elastic elongation of the suspension cords, 
and for the probable unevenness of the impact surface. 
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Effect of Error In the Time at Which Skin Is Ruptured 


Shown in figure 15 is the final impact velocity versus -u 1D , to 
which correspond the design parameters given in figures 3 through 7, for 
various errors in the time of skin rupture. Note that results are given 
only for e* = 0 and would he different if £ deviated from the design 
value. These results may, however, he used when there is a deviation 
from the design impact velocity, provided 

€ T = e u(2 + € U ) , (eg = 0) 

The reader will note that the results are given for dimensionless time 
errors. One will recall from the analysis that, when At is negative^ 
u f is just the value of u when t = t r = t^ + At, hut that when At 
is positive^ up must he calculated from equation (63) . The results_for 
negative At are indicated hy solid lines and those for positive At, 
hy dashed lines. Evidently negative errors in the time of skin rupture 
(i.e., early rupture) cause the payload to experience a higher final 
impact velocity than do positive errors. As an example of the results, 
consider the case discussed above in which u 1D = -500 ft/sec, y 2E = 0.10 
and T x = 500° R, hut in which there is no deviation from the design 
temperature or velocity. The radius for that case is 5^6 feet. Therefore 

At = 5 °° . tt/s ec At „ _At 

5.6 ft 0.0112 sec 

It can thus he found directly from figure 15(h) that the final impact 
velocity of the payload for this particular case would have the values 
in the following list for the given examples of time errors. 


t 

t, 

sec 

-Uf, 
ft/ sec 

-0.15 

-0.00168 

107 

-.10 

-.00112 

72 

-.05 

-.000560 

36 

-.01 

-.000112 

7-5 

0 

0 

0 

.01 

.000112 

6 

.05 

.000560 

29-5 

.10 

.00112 

58 

.15 

.00168 

86.5 


It is seen that the time error must he very small if the final impact 
velocity is to he small. One may note, however, that At in these results 
is proportional to the sphere radius, so that if a larger radius is used 
(and, consequently, a lower maximum acceleration) the values of At in 
the above list will he larger. 
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Demonstrated by the results of figure 15 is the fact that, if the 
sphere skin is to he ruptured to prevent rebound, it is essential to 
develop a method for rupturing the skin at the precise moment the velocity 
is zero with respect to the impact surface. 


CONCLUDING REMARKS 


The performance of the inflated sphere landing vehicle, wherein a 
payload to he landed is centrally suspended hy numerous cords, has been 
evaluated with regard to allowable mass of payload package, mass of skin 
required, mass of suspension cords required, mass of gas required, and 
radius of sphere required. The analysis suggests that the inflated 
sphere landing technique is a feasible method to use, for example, in the 
terminal phase of a moon landing which is initiated at an altitude as 
great as 25,000 to 50,000 feet from the surface. 

From the results of the analysis in this report, in which the required 
mass of suspension cords has been approximately calculated, it has become 
evident that the fraction of the total mass of the system taken up by the 
mass of cords is, in general, of significant magnitude. It is seen upon 
comparing the results with those in a previous study, in which the mass 
of cords was assumed to be negligible, that the payload mass given in that 
study closely agrees with the sum of payload and required cord mass as 
given in the present study. The results of the present report allow 
determination, however, of the actual portion of the mass available for 
payload . 

It is indicated that the skin— strength parameter and the strength- 
mass ratio of the suspension-cord material are important in determining 
the performance requirements and capabilities of the system, that is, the 
performance results depend critically on the allowable values of these 
parameters. Therefore, although the landing vehicle is shown to be 
feasible with the use of presently available materials, large gains in 
performance are to be expected as stronger, light-weight, filament-type 
materials for use in construction of flexible fabrics and cords become 
available . 

The effects of deviations from design initial temperature and impact 
velocity have been studied. The conditions have been found for which, if 
the vehicle is designed to compress to a volume somewhat greater than a 
hemisphere, a successful landing can be accomplished even if the design 
inpact velocity is exceeded and/ or the design initial temperature is not 
achieved. Critical deviations are found which, if exceeded, will result 
in the payload striking the surface with high velocity. 
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It has been found that, if the sphere skin is to he ruptured to 
prevent rehound, the effects of errors in the time at which the skin is 
ruptured may he very critical unless these time errors are exceedingly 
small, in the order of 0.1 millisecond for a 10-foot diameter sphere. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 22 , i960 
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APPENDIX 

INFLATED SPHERE LANDING VEHICLE WITH 
TAPERED PAYLOAD SUSPENSION CORDS 


In the Analysis, the payload suspension cords were assumed to he of 
constant cross-sectional area. Because the cords themselves have mass, 
then, the stress is not constant over the length of one cord at a given 
instant during the impact acceleration* The maximum allowable stress then 
occurs at only the top end of the vertical cord, and the stress in the 
remainder of the vertical cord is less than the maximum allowable at the 
instant of maximum acceleration. Hence the entire mass of cords is not 
being used to maximum advantage. The optimum use of the mass of cords 
can be realized if the stress in the entire length of the cords reaches 
the maximum allowable value during the impact. If this condition is 
imposed^ the cords will be tapered, and the appropriate cord cross- 
sectional area distribution can be calculated. While possibly not as 
practical nor as easily fabricated as the system using constant-area cords, 
the system with tapered cords is studied in this appendix to show the 
advantages in payload-landing capability. 

In addition to the symbols used in the text of the report, the 
following are applicable to the analysis in this appendix: 

A C T local area of suspension cord cross section 

A* value of A C T at t) = 0 

b parameter defined by equation (AI 3 ) 

f cp T /cp 0 f , distribution of force in the cords 

(3 2 t given by equation (A20) 

cp ! tensile force in one cord (function of 0, r|, and y (or t)) 

total vertical force in the cords at a radial distance T) in 

° < 0 < jt/ 2 

The conditions and assumptions used in the text for the case of 
constant-area cords will be retained and equations will be developed for 
the case of cords of cross-sectional area which varies so that the stress 
in the entire vertical cord has a constant value at the instant of maximum 
acceleration. Equations (l) through (5), ( 13) , and (ik) in the text 
determine the motion during impact for specified values of mp/m, m s /m, 
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n^/m , and |, provided an expression for EicvAte is known. (The procedure 
discussed in ref. 1 for finding the required value of £ to realize a 
specified y 2 may also he used.) 

Calculation of the Cord Cross-Sectional Area and 
Effects of the Cord Mass and Cord Forces 


For the purpose of determining an expression for m cv , in place of 
equation (7) one must now write for the cord mass in a radial element , 
dr), and an element of 0, d0 


where Ac 1 


Ac 1 ( T)) • 


dm,, = ^ NcpcAc’sin 9 d 9 dr] 
The total mass of cords is 


me = W cPc 



r 

A c ’dT] 


(Al) 


(A2) 


An expression for the portion of the mass of cords in motion, mc V , is then 
determined hy combining equations (Al) and (A2) and substituting the 
result for dnic into equation (10) . The expression is found, at this 
point, only in terms of the area distribution of the cord cross section: 


^v 

me 


'it p-y/cos 0 


= 1 + y , a 


A c f dr] sin 0 d G 


Ac’dT) 


(A3) 


The development which follows will lead to a determination of 
A c ’(r|). In this analysis the function cp* will he defined as the tensile 
force at the radial distance T) in the cord at angle 0. (Note the 
difference between 9' and the function 9 used in the text for cords of 
constant cross section.) Then the distribution of force in the cords, 
f , is defined by 

9« = 9*(y, 9, tj) = 9 0 '(y> Tj)f(y, 9) {kb) 

The total vertical force in the cords at a radial distance t] over 9 
from zero to rt/2 is found by substituting equations (6) and {kb) into 
the following expression: 

p9= sij 2 . 

$*(t]) = / 9*cos 9 dN c (A5) 

^0=o 


Thus 

$» = i Nc9o'(y; 


(A6) 
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where p is defined by equation (20) . The force O' (r|) is being used to 
decelerate the mass of payload, the moving mass of cords below the 
horizontal diameter, and the mass of cords above a horizontal diameter 
within a radial distance q from the center. Therefore the mass being 
decelerated by <£>' (q) is 

nip + (^cy - — ^ + — N c p c J' A c ’dq 

Because all parts of the moving portion of the system have the same 
velocity and acceleration, may be written 


du 


me 


= dt + “cv - — + 2 NcPc / V dT l 


(A7) 


From equations (a 6) and (A7) 


Tc 


t - 


N C P 




me i 

+ mey — — + — N c p c 



(A8) 


Because the cord cross-sectional area is made to vary so that the stress 
in a cord is constant throughout the length of the cord at a given time , 
one may write 


9* = crc(y> ®)A c , (ti) (A9) 

The maximum stress occurs in the vertical cord (0 = 0) at the instant 
of maximum acceleration (du/dt = ng e ) and is denoted by cr C 2 • From 
equation (A8) , then, 

(n) = n 6e(™p + ™cv 2 - + \ N cPc f Ac'dqJ (AlO) 

where mcv 2 is the value of n^y given by equation (A3) when Y - Jz- 
The solution of the integral equation, (AlO) , is 

A c *(q) = (All) 

where 

q = T]/r (A12) 


and 


and where 
to be 


= Se( nr ) 

l-* 2 ( &C 2 ./ Pc) 

A*, the value of Ac’ at q = 0, 


(A13) 


is found from equation (AlO) 
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A* 


2b 


N c Pc r 



+ m c 


f mCV2 
\ “c 



(Al4) 


Equation (All) shows that the cord cross-sectional area distribution 
should he exponential with radial distance in order to have constant 
stress in a cord. The quantity movs/W in equation (Al4) is determined 
hy the condition y = y 2 when the relationship for m^/n^ is known. 
Substitution of equation (All) into equation (A3) gives the required 
relationship: 


fcv _ l + y + 1 l(y) - l + y 
m e 2 2 L et> - 1 - 


(A15) 


where 


l(y) - f e b y/ cos ® sin G d Q 


a 


(A16) 


The transformation 


cos 



then leads to 


I a - 


■/y w2 


dw 


(ait) 

(A18) 


which has the solution 


I (y) = e'by - ye 13 - by log e y - by^ 

n=i 


(by) n - b n ~ 

n-n! 


(A19) 


The quantity m^/m^, as a function of y is now completely determined, 
by equations (A15) and (A19) , in terms of the parameter b. Thus, by use 
of equations (2), (5), and (A15) , the quantity Dv/m is completely 
determined for use in equation (3) when mp/m, m s /m, iric/m, and b are 
specified. The cord cross-sectional area distribution is given by equa- 
tion (All) where A* is given by equation (Alt) and where n^ys/mc 
may now be found from equation (A15) to be 


0 t = _ i + J2. 

P2 “c 2 


1 l(y g ) - 1 + yg 

2 L e* - 1 J 


(A20) 


The parameter b in the above development involves p 2 * As 
discussed <?n page 10 of this report, the final results can be given in 
terms of [ Ms ( cr c2 /p c ) ] and p 2 can be estimated. Also as discussed, as 
long as the skin tensile stress remains greater than zero at all points 
in the upper hemisphere, the skin above a horizontal diameter can be 
assumed to act as a rigid shell. Then the total elongation of each cord 
in 0 < Q < it/ 2 due to a small vertical deflection 5p of the payload 


package is 


SpCos 0 « 


f 9 T 

'o Ac'E 


dr! 


ffo > 0 in 0 < 0 < -? 


(A21) 


Because cr c is invariant with t), substitution of equation (A9) into 
(A21) gives the result 


Spcos 0 


a c r 

E 


a s > 0 


o < e < | 


(A22) 


and further use of equation (A9) with (A22) gives 


cp* 



(A23) 


from which equation (A4-) gives 


f = cos 6 , 

( Cg > 0 , 

o < e < jt/2) 

(A2k) 

Then, from equation (20) , if 

ff s > 0 at 

all 9 in 0 < 9 < rt/2, 



1 

q = - 


(A25) 
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As explained on page 11, for those cases in which (ct s ) 0 is reduced 
to zero at some time during the impact, the value of \i 2 will lie 
between l/3 and l/2. 


Landing Performance 


In order to calculate the design payload-landing performance, the 
required relationships among the masses of the various components of the 
system and the sphere radius are found by a procedure similar to that 
used in the text for constant-area cords. Two of the required relation- 
ships are supplied by equations (l6) and ( 38 ) in the text. A third 
expression is found by substituting equation (All) into (A2) to obtain 

me = N - c - P - crA * (e* - 1) (A26) 

b 

which may be combined with equation (Al4) to yield 


1 
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where m^ ) indicates m( )/m. A fourth relationship, obtained from 
equation (34), may be written as 

3335lflcr 

nr = — — 2 — 

mp + P 6 m s + f3 2 »mc 

The four equations, (l 6), ( 38 ), (A27) , and (A28) , may be used to find 
the four quantities m s , me, 5g, and nr in terms of nip and the various 
other input parameters involved. The product nr is implicitly involved 
in these equations since, from equation (A13) , 

b »-SE 

2Pi 

The explicit relationships can therefore be most easily found by a simple 
"trial and error" or iteration procedure . For a given value of 5L 
and given values of the various input parameters involved, an arbitrary 
value for the product nr can be used for the first trial in the itera- 
tion. The corresponding m,, can then be calculated from equation (A27) , 
and nig and 5g can be calculated from equations ( 38 ) and (l6) , respec- 
tively. Then, when the correct value of nr has been tried, equation 
(A28) will be satisfied, and the corresponding correct values of me, i5 s , 
and ihg will be found. The design quantities (i.e., the mass ratios and 
the sphere radius) then correspond to a certain impact velocity u 1; which 
can be found only by integrating the differential equation of motion, ( 3 ), 
to obtain the required g (by the procedure outlined in ref. l) and by 
making use of equation (44) . 

The equations in this analysis were programmed for solution on an 
IBM type 704 electronic data-processing machine. The differential 
equations involved were integrated by the Adams-Moulton predictor-corrector 
method. A comparison of the payload-landing capability of a sphere with 
exponentially tapered payload-suspension cords with that of a sphere 
employing constant-area cords is shown in figure l6. The results are 
shown only for y 2 = 0, cr sl /p s = 10 s ft 2 / sec 2 , and Ti = 500° R, and 
correspond to values of [P 2 (crc 2 /Pc) ] as indicated. It is seen that the 
use of the exponentially tapered cords affords a somewhat greater 
payload-landing capability at a given impact velocity. The principal 
advantage is realized when low values of mp/m are used, in which cases 
the impact velocity allowed is significantly higher. Specifically, for 
p 2 ( cr c 2 /p c ) = 1.6xl0 6 ft 2 /sec 2 (which corresponds to cr c 2 /p c = 4x10 s 
ft 2 /sec 2 if p 2 = 0.4) . a sphere with mp/m of one-fifth can land at 11 
percent higher velocity, and if mp/m is l/lO, the increase in impact 
velocity afforded is 17 percent . 


(A29) 


(A28) 
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Figure 1.- Inflated sphere during impact. 
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Figure 2.- Drop altitude above moon and corresponding impact velocity. 
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(a) Tie = 500° R; fe^ca/pc) = 1-6X10 6 ft 2 /sec 2 

Figure 3 .- Mass of payload package (hydrogen inflating gas; p a = 0; 

dsi/ps = 106 ft 2 / sec 2 ) • 
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Cb) T 1d = 200° R; p 2 (o c 2 /p c ) = 1.6x10 s ft 2 /sec 2 
Figure 3*- Continued. 
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(c) T 1D = 500° R; fe(<Wp c ) = 1-2X10 6 ft 2 / sec 2 
Figure 3-“ Continued. 
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~ u l n > ft/sec 

(<*) t 1d = 200° R; p 2 (cr c2 /p c ) = 1.2X10 6 ft 2 /sec 2 


Figure 3«- Concluded. 
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(a) T 1D = 500° R; MeUca/pc) = 1-6X10 6 ft 2 / sec 2 

Figure 4 .- Mass of skin required (hydrogen inflating gas; P a = 0; 

tfsi/Ps " 106 ft 2 / sec 2 ) • 
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(b) T 1D = 200° R; \jl 2 (o C 2 /p c ) = l-6xl0 6 ft 2 /sec 2 
Figure 4 .- Continued. 
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(c) T]_j) = 500° R; MaC^ca/Pc) = l*2xlO s ft 2 /sec 2 


Figure 4.- Continued. 
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(d) T 1D = 200° R; H 2 (crc 2 /Pc) = 1-2X10 6 ft 2 /sec 2 


Figure 4 .- Concluded. 
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(a) T 1d = 500° R; M. 2 (cr C 2 /Pc) = l«6xl0 6 ft 2 /sec 2 


Figure 5.- Mass of suspension cords required (hydrogen inflating gas; 

p a = 0; asi/ps = 106 ft 2 / sec 2 ) . 
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Cb) T 1D = 200° R; M. 2 (tfc 2 /Pc) = 1.6x10 s ft 2 / sec 2 
Figure 5 Continued. 
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(c) — 500° R; 42 ( ^02/ Pc ) = 1.2x10 s ft 2 / sec 2 


Figure 5.- Continued. 
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(d) T 1D = 200° R; ds^ca/pc) = 1-2X10 6 ft 2 /sec 2 


Figure 5 .- Concluded. 
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(a) T 1D = 500° R; Me(crc2/pc) = l-6xl0 6 ft 2 / sec 2 

Figure 6.- Mass of gas required (hydrogen inflating gas; p a = 0; 

dsi/Ps = 1° 6 ft 2 /sec 2 ). 
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0) T 1D = 200° R; |i 2 (<%>/pc) = l-6xl0 Q ft 2 / sec 2 
Figure 6.- Continued. 



(c) T 1D = 500° R; ^ 2 (^ 2 / Pc) = 1-2X10 6 ft 2 /sec 
Figi are 6.- Continued. 
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(d) T 1D = 200° R; Me(tf cz/pc) = 1.2x10 s ft 2 /sec 2 
Figure 6.- Concluded. 
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(a) Tn> « 500° R; M-2 ( ^C2 / Pc ) = 1.6x10 s ft 2 / sec 2 

Figure 7.- Radius of sphere required (hydrogen inflating gas; p a = 0; 

crsi/ps = 10 6 ft 2 /sec 2 ) . 
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-u-jjqj ft/sec 

(b) T 1D = 200° R; [i 2 (a c2 /p c ) = 1.6x10 s ft 2 / sec 2 


Figure 7*- Continued. 
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Figure J.- Continued. 
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(d) T 1D = 200° R; d 2 (ffc 2 /pc ) = 1-2X10 6 ft 2 /sec 2 
Figure 7 • - Concluded . 
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Figure 8 .- Payload at two initial temperatures (hydrogen inflating gas 
P a = 0; feCOca/pc) = 1*6x10® ft 2 /sec 2 ; a s i/ps = 10® ft 2 / sec 2 ; y 2 p = 0) 



Figure 9 .- Comparison of present analysis with preliminary analysis 
(hydrogen inflating gas; p a = 0; Me( Gq.2.1 Pc) = 1 • 6xl0 6 ft 2 /sec 2 ; 
tfsi/ps = 10 6 ft 2 /sec 2 ; Ti = 500° R; jz - 0) . 
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Figure 10.- Fractional errors in £ corresponding to fractional errors in 

ui and Tj.. 
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(a) T 1D = 500° R; p 2 (cr ca/pc) = 1.6x10 s ft 2 / sec 2 

Fig lore 11.- Final residual velocity for errors in temperature or velocity 
(hydrogen inflating gas; p a = 0; cr sl /p s = 10 s ft 2 /sec 2 ). 
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(Id) T 1d = 200° R; n 2 (o c2 /p c ) = 1.6x10 s ft 2 /sec 2 
Figure 11.- Concluded. 
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(a) T 1d = 500° R; \i 2 ( cr c2 /p c ) = 1.6x10 s ft 2 /sec 2 ; y 2 p = 0 

Figure 12.- Final location of payload package for errors in temperature or 
velocity (hydrogen inflating gas; p a = 0; cfsi/Ps = 1C|6 ft 2 / sec 2 ). 
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Figure 12.- Continued. 
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~ u i n } ft/sec 

(c) T 1D = 500° R; V-2(o c2 /p c ) = 1.6x10 s ft 2 /sec 2 ; y 2D = 0.20 


Figure 12.- Continued. 
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(d) T 1D = 200° R; m. 2 (°c 2 /pc) = 1.6x10 s ft 2 / sec 2 ; y 2 D = 0 


Figure 12.- Continued. 



(e) T lr) = 200° R; Me (o' C 2 /Pc) = l«6xl0 6 ft 2 / sec 2 ; y 2D = 0.10 


Figure 12.- Continued. 
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-u^p, ft/sec 

(f) T 1D = 200° R; |i 2 (a c2 /pc) = 1.6x10 s ft 2 /sec 2 ; y 2 D = °* 20 


Figure 12.- Concluded. 
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(a) T 1D = 500° R; H 2 (cr C 2 /pc) = 1.6xl0 6 ft 2 /sec 2 ; y 2D = 0 

Figure 13. - Final acceleration corresponding to errors in temperature or 
velocity (hydrogen inflating gas; p a = 0; Csi/ps = 10 s ft 2 /sec 2 ). 
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(t>) T 1D = 500° R; Me(<Wpc) = 1.6x10 s ft 2 /sec 2 ; y 2D = 0.10 


Figure 13*- Continued. 
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ft/sec 

(c) T]_ d = 500° R; M -2 ( ^c 2 /pc) = l*6xl0 6 ft 2 /sec 2 ; y 2 D = 0,20 


Figure 13 •- Continued. 
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(d) T 1D = 200° R; ^2(^2/ Pc) = 1.6x10 s ft 2 /sec 2 ; y 2 j) = 0 
Figure 13*- Continued. 
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-U 1D' ft / sec 

( e ) T 1D = 200° R; [x 2 (a c2 /p c ) = 1.6xl0 6 ft 2 /sec 2 ; y 2D = 0.10 


Figure 13*- Continued. 
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(f) T 1D = 200° R; h 2 (cJc 2 /pc) = l-6xl0 6 ft 2 / sec 2 ; y 2]D = 0.20 . 

Figure 13 • — Concluded • 
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(a) T 1D = 500° R; ^ 2 (^ 02 / Pc) = l>6xi0 6 ft 2 / sec 2 

Figure l4.~ Critical deviation in £ (hydrogen inflating gas; p a = 0; 

hsi/Ps = 10 s ft 2 / sec 2 ) . 
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Ob) T i;d = 200° R; Me( CT c2 /p c ) = 1.6x10 s ft 2 /sec 2 
Figure 14.- Continued- 
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(c) T 1D = 500° R; Me(cr c2 /p c ) = 1-2X10 6 ft 2 / sec 2 
Figure l4.- Continued. 
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(d) T 1D = 200° R; M-2 ( cr c2 / P c ) = 1 • 2x10 s ft 2 / sec 2 


Figure 1^.- Concluded. 
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(a) T 1d = 500° R; y 2D = 0 

Figure 15.- Final impact velocity for errors in time at which skin is 
ruptured (hydrogen inflating gas; p a = 0; cr S i/ Ps = 10 6 ft 2 / sec 2 ; 
|a 2 (tTc 2 /pc) = 1 • 6xl0 6 ft 2 / sec 2 ; = 0) . 
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Figure 15 •- Continued. 
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(c) T 1D = 500° R; y 2D = 0.20 


Figure 15 • - Continued. 
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(d) T 1D = 200° R; y 2D = 0 


Figure 15 •- Continued. 
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( e ) T 1D = 200° R; y 2D = 0.10 
Figure 15 •- Continued. 
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t id ~ 200 ° R ; y 2D -0.20 

Figure 15.- Concluded. 
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Figure 16 - Payload “ SinS 


